BCS-BEC Crossover of a Quasi-two-dimensional Fermi Gas: the Significance of 

Dressed Molecules 
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We study the crossover of a quasi-two-dimensional Fermi gas trapped in the radial plane from 
the Bardeen-Cooper-Schrieffer (BCS) regime to the Bose-Einstein condensation (BEC) regime by 
crossing a wide Feshbach resonance. We consider two effective two-dimensional Hamiltonians within 
the mean-field level, and calculate the zero temperature cloud size and number density distribution. 
For model f Hamiltonian with renormalized atom-atom interaction, we observe a constant cloud size 
for arbitrary detunings. For model 2 Hamiltonian with renormalized interactions between atoms 
and dressed molecules, the cloud size deceases from BCS to BEC side, which is consistent with the 
picture of BCS-BEC crossover. This qualitative discrepancy between the two models indicates that 
the inclusion of dressed molecules is essential for a mean-field description of quasi-two-dimensional 
Fermi systems, especially on the BEC side of the Feshbach resonance. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s 



I. INTRODUCTION 

The interest on low-dimensional Fermi systems has 
been recently reinvoked by the experimental develop- 
ments of cooling and trapping atoms in optical lat- 
tices p], [H, H[ and on atom chips Q. With the aid of 
tuning an external magnetic field through a Feshbach 
resonance, these techniques provide a fascinating pos- 
sibility of creating quasi-low-dimensional Fermi systems 
with a controllable fermion-fermion interaction. In par- 
ticular, the interaction between fcrmions can be tuned 
from a Bardeen-Cooper-Schrieffer (BCS) limit to a Bose- 
Einstein condensation (BEC) limit, such that the BCS- 
BEC crossover can be studied in quasi low dimensions. 
The BCS-BEC crossover has been extensively studied in 
three-dimensional (3D) Fermi systems, where a single- 
channel model [H and a two-channel model Q are both 
applied to give a consistent description around a wide 
Feshbach resonance. This agreement between single- 
and two-channel models is rooted from the fact that the 
closed-channel (Feshbach molecule) population is negli- 
gible near a wide resonance, so it will not cause any sig- 
nificant difference by taking the molecules into account 
(as in the two-channel model) or completely neglecting 
them (as in the single-channel model). The BCS-BEC 
crossover of a uniform two-dimensional (2D) Fermi sys- 
tem has also been considered in connection with high- 
T c superconductors pj, where an effective 2D Hamil- 
tonian with renormalized fermion-fermion interaction is 
employed. 

In this paper, we study the BCS-BEC crossover in a 
quasi-2D Fermi gas, first using an effective 2D Hamilto- 
nian with renormalized atom- atom interaction (model 1) 
[H, Q , and then a more general model with renormalized 
interaction between atoms and dressed molecules (model 
2) 10]. The dressed molecules mainly come from popula- 
tion of atoms in the excited levels along the strongly con- 
fined axial direction near a Feshbach resonance [loL [Tl| . 
When considering the effect of a weak harmonic trap 



in the two loosely confined dimensions under the local 
density approximation (LDA), we adapt the mean- field 
(MF) treatment to calculate the zero temperature cloud 
size and number density distribution in the radial plane. 
We find a significant difference between the two models. 
By using model 2, we show that the cloud size decreases 
from the limiting value of a weakly interacting Fermi gas 
as one moves from the BCS to the BEC side of the Fes- 
hbach resonance, and approaches to the limiting value of 
a weakly interacting Bose gas in the BEC limit. This be- 
havior is a signature of the BCS-BEC crossover in quasi 
two dimensions. On the contrary, model 1 fails to de- 
scribe this crossover behavior, but predicts a constant 
cloud size and identical density profile for all magnetic 
field detunings. This discrepancy implies that the MF 
results given by model 1 is unreliable, even at a quali- 
tative level. Given this qualitative discrepancy and the 
problem associated with model 1 for description of the 
two-body ground state of the system P2 , it is likely that 
the oversimplification is rooted in the model itself instead 
of the mean-field approximation. 

The quasi-2D geometry can be realized by arranging a 
one-dimensional (ID) optical lattice along the axial (z) 
direction and a weak harmonic trapping potential in the 
radial (x-y) plane, such that fermions are strongly con- 
fined along the z direction and form a series of quasi-2D 
pancake-shaped clouds Q. Each such pancake-shaped 
cloud can be considered as a quasi-2D Fermi gas when 
the axial confinement is strong enough to turn off inter- 
cloud tunneling. The strong anisotropy of trapping po- 
tentials introduces two different orders of energy scales, 
with one characterized by hw z and the other by hujj_, 
where uj z (lu±) are the trapping frequencies in the axial 
(radial) directions. The separation of these two energy 
scales (uj z ^S> uj±) allows us to first deal with the axial 
degrees of freedom and derive an effective 2D Hamilto- 
nian, and leave the radial degrees of freedom for later 
treatment. 
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II. MODEL 1 WITH RENORMALIZED 
ATOM-ATOM INTERACTION 

The effective 2D Hamiltonian for model 1 is obtained 
by assuming that the renormailied atom-atom interac- 
tion can be characterized with an effective 2D scattering 
length, with the latter derived from the exact two-body 
scattering physics [1, @. Thus, for a wide Feshbach res- 
onance where the Feshbach-molecule population is neg- 
ligible, we can write down an effective Hamiltonian only 
in terms of 2D fermionic operators at.o- and a^ r , with 
(pseudo) spin a and transverse momentum k 
The model 1 Hamiltonian thus takes the form 
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where ek — h 2 k 2 /(2m) is the 2D dispersion relation 
of fermions with mass m, [i is the chemical potential, 
and L 2 is the quantization area. The bare parameter 
V* s is connected with the physical one Vjf thro uen 

the 2D renormalization relation [V] 09 ] 1 = [Vip] 1 — 

L~ 2 (2ek + hu z ) 1 (ftio z is from the zero-point en- 
ergy), and Vfp = Vip(a s ,a z ) depends on the 3D scat- 
tering length a s and the characteristic length scale for 
axial motion a, = \/W (rnw z ) with the expression given 
in Ref. [1, 0, [lO] . Notice that the chemical potential pL 
can be a function of the radial coordinate r = (x, y) un- 
der LDA. In the following discussion, we choose hui z as 
the energy unit so that /x, V* , and £k = a 2 k 2 /2 become 
dimensionless. 

By introducing a BCS order parameter (also dimen- 
sionless in unit of hu) z ) A = {V^ / L 2 ) ^ k (ak,j.a_k,t), 
we get the zero temperature thermodynamic potential 
density 



A 4 - #k) , 



(2) 



where Ey = \f (ek — A*) 2 + A 2 is the quasi-particle exci- 
tation spectrum. The ultraviolet divergence of the sum- 
mation over k cancels with the renormalization term in 
[Vf s ] . The gap and number equations can be ob- 
tained respectively from dfl/dA 2 = and n = —dft/d/i 
(n = N/L 2 is the density of particles), leading to 
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Notice that Eq. ([3]) can be rewritten as F(a s ,a z ) = 
—fj, + y/ ' n 2 + A 2 , where the function F absorbs all the 



dependence on a s and a z . Thus, by substituting this ex- 
pression into Eq. Q , we get a closed form for the number 
equation, 



1 



irat 



F(a s ,a z ) 



(5) 



Now we take into account the harmonic trapping po- 
tential U(r) = (uj±/uj z ) 2 r 2 /(2a 2 ) in the radial plane by 
writing down the position dependent chemical poten- 
tial n(r) = fio — U(r), where hq is the chemical po- 
tential at the trap center. It can be easily shown that 
the spacial density profile is now a parabola, n(r) = 
{uj^/uj z ) 2 {R\ ¥ — r 2 )/(2na^.), with the Thomas-Fermi 
cloud size i?TF = V2hoo, z (lu z /uj±). By assigning the 
condition that the total number of particles in the trap is 
fixed by N = J n(r)d 2 r, the cloud size takes the constant 
value Rtf = Rbcs = \/2w 2 /w_L(A r ) 1 / 4 a z , which is inde- 
pendent on the 3D scattering length a s . In fact, as one 
varies the scattering length a s , the chemical potential at 
the trap center /j,q is adjusted accordingly such that the 
identical density profile is maintained. 

This result of a constant cloud size is obviously incon- 
sistent with the picture of a BCS-BEC crossover in quasi 
two dimensions. In fact, in a typical experiment with 
a z (~ //m) much greater than the interatomic interac- 
tion potential R e (~nm), the scattering of atoms in this 
quasi-2D geometry is still 3D in nature. In particular, 
fermions will form tightly bound pairs on the BEC side 
of the Feshbach resonance as they do in 3D. Thus, in the 
BEC limit when fermion pair size a pa ; r <C a z and binding 
energy \Ef,\ ^> fku z , the system essentially behaves like a 
weakly interacting gas of point-like bosons, for which one 
would expect a vanishing small cloud size in the loosely 
confined radial plane [1, [l3[ . 

The MF result of a finite cloud size in the BEC limit 
from model 1 indicates a finite interaction strength be- 
tween paired fermions, no matter how small they are in 
size. This statement can be extracted directly from the 
number equation (|5|). which can be written in the form 
/i = mra 2 — F(a s ,a z )/2. In the BEC limit, the second 
term on the right-hand side denotes one half of the bind- 
ing energy, while the first term indicates a finite inter- 
action energy per fermion pair since it is proportional to 
the number density. As a comparison, the actual equa- 
tion of state for fermion pairs one should expect must 
take the form as for a quasi-2D Bose gas in the weakly 
interacting limit Q 



[i B ~ 3na z a s 
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in which case the quasi-2D gas is treated as a 3D con- 
densate with the ground state harmonic oscillator wave 
function in the z-direction. 

The interaction strength between paired fermions can 
also be derived by writing down a Bose representation for 
this system, where the fermionic degrees of freedom are 
integrated out in the BEC limit [l4[ • This Bose represen- 
tation leads to a two-dimensional effective Hamiltonian 
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where the quartic term characterizes the bosonic interac- 
tion. Within the stationary phase approximation, the in- 
teraction strength g2 is calculated by the leading diagram 
of a four-fermion process with four external boson lines 
and four internal fermion propagators, leading to [l4j 
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where the inverse of effective interaction is connected 
with the 3D physical parameters through [Toj 
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Here, A = (—p 2 /m + \Eb\)xo(p) is the boson-fermion 
vertex, Xo(p) is the Fourier transform of the relative wave 
function Xo( r ) of two colliding fermions in the s-wave 
channel, and G (p,uj) = (iuo — p 2 /2m — \E b \/2)~ 1 is the 
free propagator for fermions. After summing over mo- 
mentum p and Matsubara frequency to, we can directly 
show that g2 indeed takes a constant value, being inde- 
pendent of the binding energy \Eb\ of paired fermions 
and hence the 3D scattering parameter a s . Thus, we 
conclude that the MF theory based on model 1 fails to 
recover the picture of a weakly interacting Bose gas of 
paired fermions in the BEC limit, and can not be directly 
applied to describe the BCS-BEC crossover in quasi two 
dimensions. 



III. MODEL 2 WITH INCLUSION OF DRESSED 
MOLECULES 

Having shown the problem associated with model 1, 
next we consider model 2 by taking into account the axi- 
ally excited states via inclusion of dressed molecules. As 
derived in Ref. fl(ij |. the effective 2D Hamiltonian takes 
the form (also in unit of hu z ), 

h 2 = ( £ k - m) a i,a a ^<? + {~2 + Xb ~ 2 ^) d q dq 
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where (c? q ) denotes the creation (annihilation) op- 
erator for dressed molecules with radial momentum q, 
and Xb, ctb, and Vb are the 2D effective bare detuning, 
atom-molecule coupling rate, and background interac- 
tion, respectively. These parameters can be related to 
the corresponding 3D parameters by matching the two- 
body physics |ldj . By introducing the order parameter 
A = a b (d ) /L + (V b /L 2 ) J2k (ak,i a -k,T)7 we obtain the 
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Here, U p = 47ra& g /a z , g 2 = n co WU p /(hw z ), and 7 P = 
fi C o(B — Bo) I (huj z ) are 3D dimensionless physical param- 
eters, where db g is the background scattering length, \i co 
is the difference in magnetic moments between the two 
channels, W is the resonance width, and B is the reso- 
nance point. The functions in Eq. (|12p take the form 



S p (x) = 

(T p [x) 
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(14) 



where r(x) is the gamma function. 

Using this model 2 Hamiltonian, we first consider a 
uniform quasi-2D Fermi gas with a fixed number density 
n, where the inhomogeneity in the radial plane is ne- 
glected. In this case, the gap and number equations JTU]) 
and (fTTj) need to be solved self-consistently for a given 
magnetic field. A typical set of results for both 6 Li and 
40 K are shown in Fig. [TJ indicates a smooth crossover 
from the BCS (right) to the BEC (left) regimes. Here, 
results obtained from model 2 (black) are compared with 
those from model 1 (gray). In this figure and the follow- 
ing calculation, we use the parameters a bg = — 1405ao, 
W = 300 G, = 2[i B for 6 Li, and a bg = 174a , 
W = 7.8 G, fi co = 1.68^s for 40 K, where ao and \xb 
are Bohr radius and Bohr magneton, respectively. 

There are two major points that need to be empha- 
sized in Fig. [T] First, when plotted as functions of the 
inverse of 3D scattering length a z /a s , the results for 6 Li 
(solid) and 40 K (dashed) are very close, manifesting the 
near resonance universal behavior. Second, the results 
from model 1 and model 2 are significantly different, es- 
pecially on the BEC side of the resonance. In particular, 
the dressed-molecule fraction in model 2 is already siz- 
able (~ 0.16) at unitarity, and becomes dominant on the 
BEC side of the resonance (see Fig.[T]:). This result pro- 
vides another signature of inadequacy of model 1 , where 
the dressed-molecule population is always assumed to be 
negligible. 
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FIG. 1: The BCS-BEC crossover behavior of a uniform quasi- 
2D Fermi gas at zero temperature, showing (a) the chemical 
potential /i, (b) the gap A, both in unit of hui z , and (c) the 
dressed-molecule fraction rib/n. Notice that the results for 6 Li 
(solid) and those for 40 K (dashed) almost coincide as plotted 
as functions of a z /a s , indicating a universal behavior around 
the resonance point. Furthermore, significant difference be- 
tween model 1 (gray) and model 2 (black) can be observed in 
(b) and (c), which shows that model 1 is oversimplified at uni- 
tarity and on the BEC side of the resonance. The parameters 
used in these plots are uj z — 2tv x 62 kHz, and na 2 z = 0.001. 



Next, we impose a radial harmonic trap U(r) and cal- 
culate the Thomas-Fermi cloud size for a fixed number 
of particles in the trap N — J 2itn(r)rdr , as shown in 
Fig. O The most important feature of Fig. [2] is that the 
cloud size given by model 2 (solid) is no longer a constant 
as predicted by model 1 (dashed). On the contrary, by 
crossing the Feshbach resonance, the cloud size decreases 
from the limiting value -Rbcs of a nonintcracting Fermi 
gas in the BCS limit, and approaches to the 3D results 
(dotted) in the BEC limit. This trend successfully recov- 
ers the corresponding physics in both the BCS and the 
BEC limits. In addition, we also find that for a given 
number of particles in the trap, the curve trend is in- 
sensitive to the radial trapping frequency u>± within the 
experimentally accessible region. (The ujj_ = 2ir x 10 
Hz and 2tt x 50 Hz results, not shown, coincide with the 
2ir x 20 Hz line and are hardly distinguishable within the 
figure resolution.) Considering the fact that there is a 
scaling relation between ui± and N such that the physics 
is only determined by N(lo±/lu z ) 2 , this insensitivity with 
respect to the radial trapping frequency suggests that 
the experimental measurement has a rather wide range 



0.2 



model 1 




-200 -100 100 

B-B„ (G) 



200 



FIG. 2: The Thomas-Fermi cloud size of a quasi-2D Fermi gas 
of 6 Li over a wide BCS-BEC crossover region. Here, results 
from model 2 (solid) are compared with those from model 1 
(dashed). All curves are normalized to the cloud size of a 
noninteracting Fermi gas -Rbcs- Notice that the results of 
model 2 recover the correct pictures in the BCS and BEC 
limits, in clear contrast to the model 1 prediction of a flat 
line. Parameters used for these two plots are u> z — 2tt x 62 
kHz, uj± = 2-7T x 20 Hz, and the total particle number N = 
10 4 . For reference, the results for an isotropic 3D Fermi gas 
with the same total particle number is also plotted (dotted), 
where a single-channel model and a two-channel model are 
both incorporated to give indistinguishable predictions. 



of tolerance on the number of atoms. 

In Fig. [3] we show the number density and the dressed- 
molecule fraction distribution along the radial direction 
for various values of a z /a s . A typical case in the BCS 
regime is shown in the top panel of Fig. [3J where the 
dressed-molecule fraction is vanishingly small, and model 
1 and model 2 predict similar cloud sizes and number 
density distributions. The middle panel shows the case 
at unitarity. As compared with model 1, notice that the 
cloud is squeezed in model 2 and the dressed-molecule 
fraction increases to a sizable value. The bottom panel 
shows a typical case in the BEC regime, where the cloud 
is squeezed further in model 2 as the dressed-molecule 
fraction becomes significant. Notice that the results of 
model 2 successfully describes the BCS-BEC crossover, 
in clear contrast to the outcome of model 1. 



IV. CONCLUSION 

In summary, we have considered in this paper the BCS- 
BEC crossover of a quasi-2D Fermi gas across a wide 
Feshbach resonance. We analyze two effective Hamilto- 
nians and compare predictions of zero temperature cloud 
size and number density distribution in the radial plane 
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within a mean-field approach and local density approxi- 
mation. Using model 1 with renormalizd atom-atom in- 
teraction, we show that the cloud size remains a con- 
stant value through the entire BCS-BEC crossover re- 
gion, which is inconsistent with the picture of a weakly 
interacting Bose gas of fermion pairs in the BEC limit. 
On the other hand, model 2 with renormalized interac- 
tion between atoms and dressed molecules predicts the 
correct trend of cloud size variation. Based on this quali- 
tative com par ison, it is likely that the inclusion of dressed 
molecules [ToL [Tlj is essential to describe the BCS-BEC 
crossover in quasi low dimensions. 
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FIG. 3: (Color online) The in-trap number density (the solid 
lines) and dressed-molecule fraction (the dashed lines) distri- 
bution along the radial direction of a quasi-2D Fermi gas of 
6 Li, obtained from model 2 (a-c) and model 1 (d-f). The top 
panels correspond to the case of a z /a s = —I (BCS side), the 
middle panels to the case of a z /a s = (unitarity), and the 
bottom panels to the case of a z /a„ = 1 (BEC side). The 
parameters are ui z = 2ir x 62 kHz, ui± — 2n x 20 Hz, and 
TV = 10 4 . 
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